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ABSTRACT 

Recent infrared observations have revealed presence of compact (radii < Rq) debris disks around 
more than a dozen of metal-rich white-dwarfs (WD), likely produced by tidal disruption of asteroids. 
Accretion of high-Z material from these disks may account for the metal contamination of these WDs. 
It was previously shown using local calculations that the Poynting-Robertson (PR) drag acting on the 
dense, optically thick disk naturally drives metal accretion onto the WD at the typical rate Mpp ss 10 8 
g s" 1 . Here we extend this local analysis by exploring global evolution of the debris disk under the 
action of the PR drag for a variety of assumptions about the disk properties. We find that massive 
disks (mass > 10 20 g), which are optically thick to incident stellar radiation inevitably give rise to 
metal accretion at rates M > Q.2Mpp. The magnitude of M and its time evolution are determined 
predominantly by the initial pattern of the radial distribution of the debris (i.e. ring-like vs. disk-like) 
but not by the total mass of the disk. The latter determines only the disk lifetime, which can be 
several Myr or longer. Evolution of an optically thick disk generically results in the development of a 
sharp outer edge of the disk. We also find that the low mass (< 10 20 g), optically thin disks exhibit 
M <C Mpr and evolve on characteristic timescale ~ 10 5 — 10 6 yr, independent of their total mass. 
Subject headings: White dwarfs — Accretion, accretion disks — Protoplanetary disks 



1. INTRODUCTION. 

Recent ground based and Spitzer observations (Zuck- 
erman & Becklin 1987; Graham et al. 1990; Farihi et 
al. 2010) have revealed the presence of infrared excesses 
in spectra of a number of metal-rich white dwarfs (WD). 
These excesses have been interpreted (Jura 2003) as aris- 
ing from reprocessing of stellar emission by a compact 
disk of refractory debris orbiting the WD. Detailed spec- 
tral modeling (Jura 2003; Farihi et al. 2010) suggests 
that these disks are geometrically thin and essentially 
flat, although in some cases there is certain evidence for 
flaring or warping (Jura et al. 2007b, 2009). Spectral 
fitting also finds the disks to be predominantly optically 
thick (Jura 2003; Farihi et al. 2010), even though this 
statement may not hold in some parts of the disk (Jura 
et al. 2007b, 2009). These properties make circum-WD 
debris disks look very similar to the ring of Saturn (Cuzzi 
2010). 

Outer radii R ou t of these disks are always found to 
be < 1 R (Jura et al. 2007; Farihi et al. 2010). This 
is close to the Roche radius Rp ~ Rq, within which 
any object of normal density (p ~ 1 g cm~ 3 ) would be 
disrupted by the tidal forces of the WD. This observa- 
tion motivated the idea (Jura 2003) that compact debris 
disks owe their origin to the tidal destruction of massive 
asteroid-like bodies, which were scattered into the low- 
periastron orbits by planets that have survived the AGB 
stage of stellar evolution. 

Inner radii of disks are usually small, Ri n ~ (15 — 
40)i?,t, where R+ is the WD radius, although the ex- 
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act number depends on the details of spectral modeling 
(Jura et al. 2007). In some cases Ri n is close to the dis- 
tance at which a refractory body heated by the direct 
starlight would have a temperature equal to the subli- 
mation temperature T s , suggesting that the inner cutoff 
of the disk is caused by particle sublimation. However, 
in some cases such interpretation implies very high value 
of T s (Jura et al. 2007; Brinkworth et al. 2009; Dufour et 
al. 2010). Nevertheless, it is always safe to assume that 
the ratio of the outer to inner disk radii is rather modest, 
Rout I R%n ^ 5. This will turn out being important later 
on, see ^44.14.1.11 

A significant reservoir of refractory material in the im- 
mediate vicinity of the WD can naturally be responsi- 
ble for contaminating the atmospheres of these WDs by 
metals (Jura 2003). This idea is strongly supported by 
the fact that the infrared excesses due to compact de- 
bris disks have been observed only around metal-rich 
WDs with rather high metal accretion rates (Farihi et 
al. 2009, 2010). Alternative model of the WD metal pol- 
lution by interstellar accretion (Dupuis et al. 1993a) runs 
into problem trying to account for the abundance ratios 
of different elements in the WD atmospheres, in partic- 
ular, greatly suppressed hydrogen abundance (Dupuis et 
al. 1993b). 

The new model of metal accretion from a circumstel- 
lar reservoir of high-Z material raises a natural question: 
how do metals get transported to the WD surface from 
the debris disk, which has an inner gap with the radius 
It is reasonable to assume that at Ri n disk parti- 
cles get converted into metal gas by the intense heating 
due to the WD radiation, which then viscously accretes 
onto the WD. However, the issue of whether the debris 
disk can sustain the high metal accretion rates Mz up 
to several xlO 10 g s _1 inferred in some systems has not 
been clear for a long time. 
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Recently Rafikov (2011a; hereafter Rll) showed that 
coupling of the disk to stellar radiation via the Poynting- 
Robertson (PR) effect (Burns et al. 1979) can result in 
metal accretion rates Mz ~ 10 8 g s^ 1 . Later Rafikov 
(2011b) has also demonstrated that under certain cir- 
cumstances even higher values of Mz can arise from the 
interaction between the debris disk and metal gas that is 
produced by the sublimation at Ri n . 

The model of radiatively driven accretion of Rll did 
not cover the large scale evolution of the debris disk re- 
sulting from the PR drag, and made certain assumptions 
(e.g. high optical depths of the disk) which were not rig- 
orously verified. The goal of this work is to extend the 
analysis of Rll and to develop a detailed global model 
of the compact debris disk evolution caused by the PR 
drag. 

The paper is organized as follows. In <2]wc outline the 
basic picture of the PR-driven debris accretion and derive 
master equation (|17[) that describes global evolution of 
the disk. We then explore in both analytically and 
numerically the evolution of a low mass, optically thin 
disk of debris. In £|4]we study global evolution of massive, 
optically thick debris disks starting with different initial 
spatial distributions of debris around the WD (ring-like, 
34 4.14.1.11 T4 4.14.1.21 or disk-like, 344.14.1.3j) to see the 
effect on the global disk evolution. We discuss our results 
and their observational implications in 

2. DESCRIPTION OF THE MODEL. 

In the following we consider an axisymmetric disk of 
particles extending from Ri n to R ou t in radius. The inner 
radius may coincide with the sublimation radius R s , at 
which the effective temperature of particles equals the 
sublimation temperature T s : 
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where R+ is the WD radius, T s rj 1500 K for silicate 
grains, and T*^ = ^/(lO 4 K) is the normalized stellar 
temperature T*. Taking i?* w O.OIRq one finds R s w 
0.2 Rq, in agreement with observationally inferred inner 
radii of compact debris disks (Jura et al. 2007, 2009a). 

When particles reach the sublimation radius they pro- 
duce metallic gas, which joins the gaseous disk extend- 
ing down to the WD surface. Metal accretion onto the 
WD surface proceeds through this disk. Although this 
metal gas also spreads outward from the sublimation ra- 
dius (Melis et al. 2010) and under certain circumstances 
can substantially affect debris disk evolution (Rafikov 
2011b), in this work, following Rll, we concentrate only 
on effects associated with the PR drag. Thus, here we ne- 
glect presence of the gas exterior of R s and its interaction 
with the debris disk. According to Rafikov (2011b) this 
is a valid approximation as long as the viscous timescale 
in the gaseous disk is shorter than the time on which 
this disk can be replenished by the sublimation of de- 
bris. Then the gas does not accumulate at r « Ri„ and 
its density is always low enough for the gas drag to not 
affect the debris disk appreciably. 

We characterize the disk at each point by its surface 
density E(r) and optical depth 
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where p is the bulk density of particles and a is their char- 
acteristic size. The particle size a is not well constrained 
by the existing observations, but its actual value becomes 
important only for low mass disks with small r, see $3] 

Following Friedjung (1985) and Rll we represent irra- 
diation of the disk by the WD using a single incidence 
angle a at each radius r from the disk (the so-called 
"lamp-post" illumination model described in Rll): 
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where i?* is the radius of the star. We can then introduce 
optical depth to incoming starlight tii according to the 
following definition: 

T|| = ol~ 1 t. (4) 

This variable is very important for our subsequent analy- 
sis. In the following we will call the disk optically thick or 
thin based on whether tii and not r is greater or smaller 
than unity. 

2.1. Accretion through the disk. 

Rll demonstrated that the mass accretion rate M 
through the disk of solids driven by the PR drag may 
be written as 



M(r) = ai(rU r %. 



(5) 



Here is the WD luminosity, c is the speed of light and 
function 



V = 1 - e~ T » 



(6) 



gives the fraction of incoming starlight that is absorbed 
by the disk. 

It is instructive to examine different limits of the ex- 
pression ([5]). If the disk is optically thick to incident 
stellar radiation, tii 1 (which cannot be the case ev- 
erywhere at all times, as we will see below), then the 
accretion rate is independent of tii and is given simply 
by 



M T|| »i(r) = M(r,Tn > 1) = a(r) 



(7) 



In particular, at the sublimation radius given by equation 
([1]) this accretion rate is equal to 



M 00 =M Tll>>1 (R s ) = ^cr(^^ 
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T500K, 

Rll has shown this numerical estimate to be consistent 
with the lower envelope of the Mz values inferred for 
the metal-rich WDs exhibiting detectable IR emission 
associated with their debris disks. 

In the opposite limit of a disk, which is optically thin 
to incident stellar radiation, tii <C 1, one finds 

h 

c 2 

It is easy to show that this expression coincides with 
M that one would calculate by simply assuming disk par- 
ticles to be directly illuminated by starlight (i.e. unob- 
scured by other particles) and considering each of them 
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as independent from others. Azimuthal PR drag force 
acting on a single, perfectly absorbing particle of mass 
m = (47r/3)/oa 3 is given by 



v Anr 2 



(11) 



where Qk is Keplerian angular velocity. This force gives 
rise to radial drift velocity 
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and results in mass accretion rate M = 27rru r E, which 
is easily shown to reduce to equation (ITU1) . 

2.2. Evolution equations 

Evolution of the debris disk is described by the conti- 
nuity equation 

<9£ 1 dM 
dt 2wr dr 
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with M given by equation ([5]). 

We introduce new dimensionless time and space vari- 
ables 



where 
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is the characteristic timescale of the problem. If Ri n 
coincides with the sublimation radius R s defined in (fTj) 
then 
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where we took p = 3 g cm -3 . 

Using equations (0-© and definitions (fT4")l . (| 15[) we 
can bring the continuity equation (|13|) to a scale-free 
form: 



d_ 

dx 
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A notable fact about equation (|17p is that it does not 
contain any free parameters. 

Once T|| (x, T) is obtained by solving this equation, one 
can easily infer mass accretion rate at any point in the 
disk since 



M(x,T) = M 



1 



(18) 



where is defined by equation (|S|). Because the vis- 
cous timescale of the gas produced by debris sublimation 
is very short compared to the disk evolution time (Rll), 
the metal accretion rate onto the WD surface Mz is given 
simply by 



M Z (T) = M{x = 1,T) = M 



1 



,-T||(x=l,T) 



.(19) 



In the following we will use M Z (T) and M(x = 1,T) 
interchangingly. 



Analysis of equation (|17[) may be simplified if we 
introduce a new function y(x,T) = M(x,T)/ = 
(1 — e~ T ii)/x, which is just the mass accretion rate 

through the debris disk normalized by Mqq. Then equa- 
tion (|17l) transforms to 



dy 



(20) 
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It has an implicit solution 

y = f(xy - y 2 T + M 1 - x v)), 

where / is an arbitrary function, which is set by the ini- 
tial surface density distribution of the debris. Unfortu- 
nately, inferring / from y{x, T = 0) is not a trivial task, 
which makes analysis of solution (|21[) very difficult. Nev- 
ertheless, we can still deduce useful information about 
the behavior of this solution in some specific limits, see 
Sj3]and Appendix |A"1 

We now proceed to investigate the details of the global 
disk evolution with different initial conditions. 

3. LOW MASS DISKS. 

We start by considering the case of a tenuous disk, 
in which the optical depth for incident stellar radia- 
tion T|| <§; 1. In that case according to equation (fl"8]) 
y ~ t~\\/x <C x and the nonlinear term ydy/dx can be 
neglected in equation (|20]l reducing it to 

dy 1 dy 

df ~ xdx 

This linear equation has an explicit solution 
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(23) 



where yo{x) = y(x, T = 0) ~ T|| (x, T — 0)/x is the initial 
spatial distribution of variable y set by the initial surface 
density distribution T||(x,T = 0). 

If initially the debris was concentrated in a structure (a 
disk or a ring) with characteristic dimensionless scale x 
(physical scale Ri n x), then it follows from equation (|22|) 
that the characteristic time Uhin on which rn evolves is 
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This time is independent of the disk mass but is sensitive 
to its size x and particle properties, see equation (|16[) . 
This is not surprising since when rn 1 particles interact 
with the WD radiation independently of each other and 
drift towards the WD on a PR timescale of an individual 
particle. The latter coincides with tthin up to constant 
factors and is indeed a function of /?, a and the distance 
r that needs to be traveled. 

Disks which are optically thin to incident stellar radia- 
tion must satisfy a certain mass constraint. The mass of 
a disk extending out to R ou t = x ou tRin can be written 
as 



Rout Xoui 

M d = 2ir J E(r)rdr = M J 
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where the characteristic disk mass Mn is defined as 
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Fig. 1. — Evolution of a low mass debris disk around a WD 
with initial ry <C 1 everywhere. Initial density profile is given by 
equation l|27|l with rii rj = 0.1, xq = 5, a = 0.7. Thick curves are 
the numerical results for the optical depth to incoming starlight 
T\i(x,T) (related to E(r, t) via equations (2)-(4)) at different mo- 
ments of time indicated on the panel, while the thin curves repre- 
sent analytical solution J28j. 



and the numerical estimate assumes p — 3 g cm" 3 , and 
R± m O.Oli?©. This numerical value of Mq corresponds 
to the disruption of an asteroid with the diameter of 40 
km. 

The disk which is optically thin everywhere (tii < 1) 
should have mass Md < Mq since x ou t — 1 is of order 
unity for WD debris disks. 

3.1. Numerical results. 
We numerically integrated equation (|17p for an initial 



ring-like density profile given by 
T||(x,0) « xy (x) = T|| !0 exp 



(x - x ) 2 



(27) 



with T|| Q <C 1. In Figure Q] we show the results obtained 
for the case tii o = 0.1, xq = 5, er = 0.7, i.e. a rather 
narrow, tenuous ring, with the radius large compared to 
Ri n . We also plot analytical solution for tiiq <C 1 
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,(28) 



which follows from (f23|) . One can see that analytical solu- 
tion reproduces the numerical results quite well. As time 
goes by the ring of material drifts towards the WD and 
broadens, simply because the speed of inward migration 
due to the PR drag increases as x decreases. This even- 
tually results in substantial surface density reaching the 
inner radius and giving rise to sublimation there. Metal 
gas produced at i?i„ then directly accretes onto the WD. 

Metal accretion rate Mz <C in the optically thin 
case, simply because the initial value of tii <C 1 every- 
where in the disk and equation (|23[) then predicts that 



y(x = 1,T) -C 1 for any T. Using solution (l23l) we can 
derive an explicit formula for M(x = 1,T) for the initial 
debris density distribution (|2"T|) : 



M(x = 1,T) 



: exp 



(VT+2T-x ) 2 



(29) 



V1 + 2T 

In Figure [2] we compare this expression with the nu- 
merically derived M(x — 1,T) and again the agree- 
ment is reasonably good. Mass accretion rate peaks at 
T rj xl/2 « 12 when the expression inside the exponen- 
tial in (|29l) goes to zero. The maximum accretion rate 
achieved at this moment is about MooTiio/xo ~ 0.02A/OO, 
in agreement the numerical results shown in Figure [2] 
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Fig. 2. — Metal accretion rate at r = Ri n for a low-mass opti- 
cally thin disk (same as in Figure [Til, normalized by Moo- Solid 
and dotted curves are the numerical and analytical [equation H29H 1 
results. 



4. HIGH MASS DISKS. 

Now we consider evolution of massive debris disks with 
mass exceeding Mo defined by (f26|) . Such disks have 
T|| ^> 1 at least at some radii. This condition cannot be 
true throughout the whole disk as there should always be 
regions in which tii < 1, e.g. near the disk edges. Even 
if they do not exist initially such optically thin regions 
naturally appear as a result of disk evolution as we show 
below. 

Whenever tii 3> 1 one can neglect the term propor- 
tional to e~ T n in equation (| 1 T|) and obtain the following 
simple solution: 



T\\(x,T) = T 
or, in dimensional units, 



(x,T 



T,(r,t) = Y,(r,t 



T 



2 i?* 



(30) 



(31) 



This solution implies that in parts of the disk, which are 
optically thick to incoming stellar radiation the surface 
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density steadily decreases in time at a constant rate set 
only by the distance from the WD. This behavior is a di- 
rect consequence of the M saturation at the level M T|| >i 
independent of t\\ in the case t\\ 3> 1, see equation ([7]). 

At a glance this kind of evolution looks very strange as 
the debris does not seem to move anywhere — E simply 
goes down at all radii. This interesting behavior can 
be understood only by accounting for the optically thin 
regions that should inevitably be present in the disk and 
their interplay with the optically thick part of the disk. 
We demonstrate this in £144.11 

Solution (|30l) also implies that the characteristic 
timescale of the disk evolution is now 

t T\\{x)x 2 = — r V '- (32) 
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or TtMck — Tn {x)x 2 in dimensionless units. Clearly, tthick 
is (up to factors of order unity) just the time, on which 
the characteristic disk mass ~ E(r)r 2 gets exhausted by 
accretion at the rate M T||3> i(r). Unlike t t hm defined by 
equation (|2"4")l tthick is independent of particle properties 
— p and a, but scales linearly with the disk mass. 



4.1. Numerical 



dts. 



We numerically integrated evolution of several initially 
optically thick density distributions. The morphologies 
we consider range from a narrow ring to an extended disk 
of debris. They represent different initial spatial distri- 
butions of the debris that might possibly result from the 
tidal disruption of an asteroid. We choose all disks to 
have the same total mass Md — 350Afo initially and con- 
centrate on exploring the differences of their evolutionary 
routes caused purely by the morphology. 

4.1.1. Large and narrow ring. 

We start by exploring the initial distribution in the 
form (1271) with xq — 5, a — 0.7 but now we take Tiirj = 
280 so that the disk mass M d = 350M > M . Resulting 
evolution of tii is shown in Figure [3] Apart from the 
theoretically expected steady decrease of tii with time 
in the bulk of the disk where t\\ 3> 1 one immediately 
notices two interesting features of the solutions: (1) a 
sharp drop of the disk surface density at the outer edge of 
the disk, and (2) an optically thin (t\\ < 1) tail of debris 
extending from the optically thick region all the way to 
Ri n , which remains stable throughout the evolution of 
the disk. We now discuss the nature of these features in 
more detail. 

Sharp outer edge appears because the evolution 
timescale of a particular disk region tthick oc tii , see equa- 
tion (|32|) . Initial density distribution in the form (|27|) has 
T|| increasing inwards in the outer part of the ring. This 
causes the outermost extremity of the ring (which is op- 
tically thin and evolves on timescale Uhin) to migrate 
towards the WD faster than the parts of the ring closer 
to the WD do. As a result, the outermost disk material 
catches up with the debris at smaller x having higher tii , 
leading to a debris pileup there and formation of a sharp 
outer edge. This edge appears on a timescale ~ Uhin 
since this is the characteristic time, on which the ma- 
terial from the optically thin outer disk region reaches 
x xq where most of the ring mass is concentrated. We 
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Fig. 3. — Evolution of tii for an optically thick disk with the 
gaussian initial density profile H27I I. assuming xq = 5, a = 0.7 and 
total mass = 350Mo (large and narrow ring) . Radial profiles of 
T|| are shown at different moments of time indicated on the panel. 
Note the appearance of an optically thin tail of debris extending 
to x = 1 and the sharpness of the outer disk edge. 



provide a more rigorous justification for the outer edge 
appearance in Appendix [X] 

Origin of the quasi-stationary tail with T|| < 1 at small 
r can be understood as follows. Even if initially the disk 
has T|| 3> 1 everywhere such tail would rapidly form since, 
according to equation f[30]) , at any distance from the WD 
T|| should switch into the optically thin regime within fi- 
nite time. This happens fastest at the inner edge of the 
disk. After the tail has formed it evolves on a timescale 
tthin given by equation (|24[) . Inner parts of the tail evolve 
more rapidly than the outer ones and a quasi-steady state 
is attained in the tail with M almost independent of r. 
This point is illustrated in Figure HI where we show the 
run of M{x)/Moa for a calculation displayed in Figure[3] 
at different moments of time. One can easily see that af- 
ter some evolution has taken place M(x) indeed becomes 
independent of x in the tail. 

Thus, the behavior of tii (x, T) in the tail can be de- 
scribed by the equation (cf. equation ([T51l ) 



1 - exp [-t\\(x,T)] 



= DtailiT) = 



M ta u(T) 

Moo 



(33) 



where Mtail (t) is the mass accretion rate through the tail 
(and ytaii is the dimensionless analog of this quantity), 
which in general is a function of time. Equation (|33p 
describes the behavior of tii in Figure [3] as a function of 
x quite well. 

At some radius x T= \ = r T= \/ Ri n the tail becomes op- 
tically thick, i.e. T||(x T= i) = 1. A this point the tail 
merges with the optically thick part of the disk and the 
evolution timescale, which is now given by tthicki rapidly 
increases, see equation (|3"2")l . Accretion rate through the 
tail M t aii{T) is set almost exclusively by the value of 
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Fig. 4. — Radial distribution of the dimensionless mass accretion 
rate y(x,T) = M(x,T)/Moa (equation (15) 1 through the debris 
disk at different moments of time. Results are shown for the same 
disk as in Fig. \3\ Note that y — ¥ const near x = 1, i.e. in the 
optically thin tail. 



x T —i since by definition 1 — exp [— rn (x T —i, T)l 
dependent of T, so that 
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x T =i(T) 
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(34) 



at all times. Thus, it is the dependence of x T —i on T that 
determines how the mass accretion rate through the tail 
and at Ri n changes with time. 

To find x T= i(T) we need to set the left hand side of 
equation to unity and solve the resultant equation 
for a: as a function of T. Since here we consider the case 
of T\\{x,T = 0) 3> 1 the relation between x T —\ and T is 
provided with sufficient accuracy by solving the equation 



n(x T =i,T = 0) 



(35) 



for a given initial density distribution in the disk. 

For the Gaussian initial distribution in the form ([27 
one finds that 
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(36) 



Whenever xq ^> 1 and a <C £o, which is a reasonable 
approximation for the case shown in Figure [3l equation 
(|3"6"|) predicts that x T= \ « xo and that x T —\ depends on 
time T only weakly (logarithmically). This explains the 
time invariance of the optically thin tail in Figure [3] and 
only weak time dependence of M in the tail in Figure 01 
the mass flux through the tail 

Mum « MooXq 1 (37) 

is almost independent of T if the initial mass distribution 
has the form of a large and narrow ring (i.e. a <C xq). 



As a result, the flux of metals onto the WD M(x = 
1, T) also varies with time only weakly for ring-like initial 
distribution of the debris, as can be seen in Figure[5] For 
debris disks around WDs produced by tidal disruption 
of asteroids xq <"o and so equation (|37|) guarantees that 
Mz stays above 0.2Moo or so, see Figured] 
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Fig. 5. — Mass accretion rate at the inner disk radius M(x = 
1,T) normalized by Moo for different ini tial distrib utions of the 
debris surface density: large narrow ring f i|4 4,l 4.1.11 dashed line), 
s mall ring (j!4 4.1 4.1.21 dot-dashed line), and extended top- hat disk 
( !|4 1.1 4.1.31 solid line). Dotted curve shows analytical prediction 
H39H for Mz in the top-hat case. Disk mass is the same in all three 
cases: M d = 350M . 



Equation (|33|) also predicts that rn does not go signif- 
icantly below unity even at the inner edge of the disk: 
n(x — l,T) ~ Xq 1 depends only on xq, see equation 
(I33p . and does not change appreciably with time for a 
narrow ring. 

It is clear that the general evolutionary picture pre- 
sented in Figure [3] would still hold if we assume a ring 
density profile other than Gaussian, as long as the width 
of the ring is small and its radius is large compared to 

4.1.2. Small ring. 

We now consider evolution of initial distribution in the 
Gaussian form (|27[) but with xq = 2, a = 0.7, and trun- 
cated at x = 1 (i.e. Tii = for x < 1), while keeping 
the total disk mass the same, Md — 350M . In this case 
the width of the ring is no longer small compared to its 
radius. Snapshots of the tii profile at different time T are 
presented in Figure [HI 

Since xq — 1 ~ a this distribution has significant rii ~ 
40 at x = 1 at the very beginning. According to equation 
((33)) it takes time T r\\ t0 (x = 1) « 40 for t\\(x = 1) 
to drop to small value ~ 1, and during this initial stage 
M(x — 1,T) sa Mqo, as Figure [5] clearly shows. 

Beyond this point a marginally optically thin tail de- 
velops near x = 1 but it never reaches a steady state. 
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This can be understood by inspection of the transcenden- 
tal equation (|36|) determining the crossover radius x T =i 
(and the M through the tail): whenever the width of 
the ring a is not small compared to its radius xq this 
equation does not have a time-independent solution. As 
a result, Mz = M{x = 1,T) noticeably varies in time, 
see Figure [5] 

Also note that the sharp outer edge of the disk is again 
present in Figure [5] — its origin is the same as in the case 
of the large and narrow ring explored in £144.1 4.1.11 



1000 



E- 



10 



0.1 



: 1 




I I 
T = 


: 


- //ft 1 

://// i 

- ;/' 
" M 


vv 

V 


T = 

\ T = 
\ T = 


ioo ; 

200 
300 i 


-> 

;; 

■'ji i 








¥ ! 








r 1 








; i 

: i 
I 

■ 1 




i \ 





3 

X 



Fi g. 6. — Same as Figure [3] but for the gaussian initial profile 
l|27|l with xo = 2, a = 0.7 and total disk mass = 350Afo (small 
ring). 



4.1.3. Extended disk. 

Finally we explore the case of an extended initial den- 
sity distribution in the form of a disk described by a 
"top-hat" initial condition 

T-jj (x, T = 0) = T|| )0 = 100, xi < x < x 2 , (38) 

where x\ = 1.5 and xi — 5. Evolution of this distribution 
is shown in Figure [7] and one can again clearly see the 
development of an optically thin tail at the inner edge of 
the disk and the sharp outer edge. 

The tail appears on a timescale T t hin{x\) ~ x\ ps 2. 
Initially this tail attaches directly to x T= \ — X\ because 
it takes time T^^ixi) ~ T \\,o x \ ~ 220 (see equation 
(|3"5")) ') for T||(xi) to go down to unity. During this pe- 
riod of time Mz ~ %i Moo ~ 0.66-Mqo stays virtually 
constant. These simple analytical conclusions are in ex- 
cellent agreement with the numerical results presented in 
Figure 

After rii(xi) has dropped below unity the crossover ra- 
dius x T= i starts moving out and M tai i goes down. Using 
equations (|3"5j) and (|3"8"]) we infer that in the top-hat case 
x T =i w (T/r^o) 1 / 2 for T > T t hick{x\), resulting in 

1/2 



M z w M 



(?) 



This analytical prediction is shown by dotted line in Fig- 
ure [5] and clearly agrees well with the numerical result 
for the initial top-hat density distribution. 



1000 




Fig. 7. — Same as Figure [3] but for a "top-hat" initial density 
distribution rn(x,T = 0) = 100 for 1.5 < x < 5. Note that both 
the optically thin tail near x = 1 and the sharp outer edge are still 
present. 



5. DISCUSSION. 

In this paper we focused on studying how the debris 
disk evolves due to the PR drag alone. Among other 
radiative effects potentially affecting such disks we can 
mention the direct transfer of the radiative angular mo- 
mentum to the disk that occurs if the WD is rapidly 
rotating. Previously Walker & Meszaros (1989) sug- 
gested this process to be important for the evolution of 
accretion disks around neutron stars experiencing Type 
I X-ray bursts, which often rotate quite rapidly. In the 
case of WDs there is no good evidence for rapid rotation 
(Kawaler 2004) with several disk-hosting systems show- 
ing v ro tsmi < 20 km s _1 (Gansicke et al. 2006, 2007). 
However, one can easily demonstrate that even if WD 
were to rotate at breakup speed the radiative angular 
momentum transport would still not produce WD disk 
evolution comparable to that caused by the PR drag. 

Indeed, the angular momentum carried by photons 
emitted by the WD rotating with equatorial speed v ro t 
can be estimated as H ~ (L ir /c 2 )R ir v rot . Disk absorbs 
some fraction of this angular momentum which exerts ra- 
diative torque T ra d on the disk. The ratio of this torque 
to the torque due to the PR drag Tpr is easily shown to 
be 



rail 



R*v r 



T P 



(40) 



(39) 



'PR SIkt 2 ' 

where Qk is the Keplerian angular frequency. Even if the 
WD rotates at breakup and v rot ~ Hk(R*)R* one still 
finds Trad/Tpp. ~ (Rt/r) 1 ' 2 , which is less than unity be- 
cause the debris disk lies quite far from the WD surface, 
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i.e. r>fi i . Moreover, measured rotational speeds of the 
aforementioned WDs are ~ 10 -2 of their breakup speed, 
so the direct radiative transfer of angular momentum to 
debris has negligible effect on disks in these systems. 

One of the key results of this work is that the low metal 
accretion rates onto the WD Mz -C Moo ar e possible 
only if the debris disk hosted by the WD has small mass, 
below Mo given by equation (|26|) . and is everywhere op- 
tically thin to incident stellar radiation. This is demon- 
strated in $3] where we find Mz ~ MooT||(i = 1) <C Mx, 
(see equation (|23p ) throughout the whole evolution of the 
low mass disk. 

On the contrary, if the disk is massive enough to 
contain optically thick regions, which typically requires 
M d > 10 20 - 10 21 g, then M z is guaranteed to not fall 
below (Ri n / Rout) Mqc ~ 0.2Moo, as we have shown in <21 
This is in agreement with the statement in Rll that Mz 
should not deviate significantly from M^ for optically 
thick disks. 

This lower bound on Mz for massive disks is very ro- 
bust, in particular it does not depend on the disk mass. 
This point is illustrated in Figure [5] where we show the 
time evolution of Mz for three disks with the same initial 
Gaussian density distribution (|27j) but with different ini- 
tial masses Md set to be in the ratio 0.5 : 1 : 2. One can 
easily see that the only thing, which is different between 
these three cases is the disk lifetime (which is always pro- 
portional to the disk mass), while Mz is essentially the 
same. 
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Fig. 8. — Time evolution of the mass accretion rate at the inner 
radius M(x = 1, T) for three massive, opt ically thick disks with the 
same initial gaussian density profile H27I I and xq = 5, cr = 0.7, but 
different initial masses in the ratio 0.5:1:2. Note that the amplitude 
of M(x = 1, T) is essentially independent of the disk mass, only 
the disk lifetime depends on it. 

The time evolution of Mz in the optically thick case is 
found to depend predominantly on the shape of the initial 
density distribution. Debris distributions in the form of 
a narrow ring tend to produce Mz only weakly varying 



with time. Extended initial distributions of the debris 
with radial width comparable to mean radius generally 
lead to Mz varying in time and decaying by a factor of 
several by the end of the disk lifetime, see Figure [5j 

Based on these results one should expect WDs with 
strong IR excesses, signifying massive, optically thick de- 
bris disk around them, to exhibit Mz > 0.2Moo. This is 
exactly what the current data seem to suggest: no WD 
with a debris disk detected via its IR signature is inferred 
to have M z below M x (Rll). 

However, if the systems exhibiting IR excesses were 
found to show Mz <C , the debris disks around them 
should be optically thin. It is then natural to expect a 
positive correlation between the strength of the IR ex- 
cess and the Mz in such systems, since both are linearly 
proportional to rii <C 1 in the optically thin regime. In- 
deed, Mz oc T|| according to (f2"3")l while the amount of the 
WD luminosity that the optically thin disk intercepts and 
reradiates in the near-IR is 

L^^f^-r^dr, (41) 

which also scales positively with rii . 

The timescale on which debris disks get depleted can 
be quite different in the optically thick and thin cases. 
Characteristic lifetime of an optically thin disk composed 
of ~ 1 cm particles t t hm is several 10 yr, as equations 
(fl6|) and ([24]) demonstrate. This timescale depends only 
on the debris particle properties — size and density, be- 
ing longer for bigger and denser particles — and the outer 
extent of the disk (for R ou t/Rin ~ 5 one finds disk life- 
time approaching a Myr). 

In the case of an optically thick disk the lifetime can 
be considerably longer since tthick scales linearly with the 
disk mass, see equation (|32[) . Thus, massive disks require 
long time to get exhausted by the PR drag alone. For 
example, an optically thick disk created by the tidal dis- 
ruption of an asteroid with the diameter of 300 km would 
have a mass of 4 x 10 22 g (assuming bulk density of 3 g 
cm~ 3 ). Accretion at the steady rate Mz = 3 x 10 7 g s _1 
(pa 0.5Moo given by the numerical estimate in equation 
((H)) would deplete this disk only within 40 Myr. Such 
massive disks, however, are likely to be affected by the 
interaction with the gas that is produced by sublimation 
of solid debris at the inner edge of the disk. They may 
then evolve in a runaway fashion as described in Rafikov 
(2011b) and get exhausted much earlier. 

An interesting finding of this work is the universal ap- 
pearance of the sharp outer edge in optically thick disks, 
caused by the faster inward migration of the disk annuli 
having smaller optical depth. This outer edge steadily 
moves inwards which, combined with the outward ex- 
pansion of the optically thin tail at small radii, leads to 
the gradual reduction of the radial width of the optically 
thick part of the debris disk. Thus, even an extended 
disk can turn into a narrow ring with time as can be 
easily seen in Figure 

Observational evidence for sharp outer edges would 
strongly support the picture of disk evolution due to the 
PR drag outlined in this work. Interestingly, Jura et al. 
(2007b, 2009) have previously claimed the need for ex- 
istence of the outer optically thin regions in debris disks 
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around some WDs. However, inferring the sharpness of 
the outer disk edge and the precise value of its optical 
depth is a nontrivial task since any knowledge about the 
radial density distribution in these disks is based on mod- 
eling their spectral energy distributions. Inverting disk 
spectrum to determine the fine details of the radial dis- 
tribution of debris is a rather ambiguous procedure, ad- 
ditionally complicated by the possibility of disk warping 



or flaring (Jura et al. 2007b, 2009). 



The financial support for this work is provided by the 
Sloan Foundation and NASA via grant NNX08AH87G. 
KVB and RRR thank Princeton University and Lebedev 
Physical Institute for hospitality. 



REFERENCES 



Brinkworth, C. S., Gansicke, B. T., Marsh, T. R., Hoard, D. W., 

k Tappert, C. 2009, ApJ, 696, 1402 
Burns, J. A., Lamy, P. L., k Soter, S. 1979, Icarus, 40, 1 
Cuzzi, J. N. et al. 2010, Science, 327, 1470 

Dufour, P., Kilic, M., Fontaine, G., Bergeron, P., Lachapelle, F.-R., 
Kleinman, S. J., k Leggett, S. K. 2010, ApJ, 719, 803 

Dupuis, J., Fontaine, G., Pelletier, C, k Wesemael, F. 1993a, 
ApJS, 84, 73 

Dupuis, J., Fontaine, G., k Wesemael, F. 1993b, ApJS, 87, 345 
Ehrenreich, D. et al. 2011, A&A, 525, 85 

Farihi, J., Zuckerman, B., k Becklin, E. E. 2008, ApJ, 674, 431 
Farihi, J., Jura, M., k Zuckerman, B. 2009, ApJ, 694, 805 
Farihi, J., Jura, M., Lee, J.-E., k Zuckerman, B. 2010, ApJ, 714, 
1386 

Friedjung, M. 1985, A&A, 146, 366 

Gansicke, B. T., Marsh, T. R., Southworth, J., & Rebassa- 

Mansergas, A. 2006, Science, 314, 1908 
Gansicke, B. T., Marsh, T. R., & Southworth, J. 2007, MNRAS, 

380, L35 

Graham, J. R., Matthews, K., Neugebauer, G, & Soifer, B. T. 
1990, ApJ, 357, 216 



Jura, M. 2003, ApJL, 584, L91 

Jura, M., Farihi, J., & Zuckerman, B. 2007a, ApJ, 663, 1285 
Jura, M., Farihi, J., Zuckerman, B., & Becklin, E. E. 2007b, AJ, 
133, 1927 

Jura, M., Farihi, J., & Zuckerman, B. 2009a, AJ, 137, 3191 
Jura, M., Muno, M. P., Farihi, J., k Zuckerman, B. 2009b, ApJ, 
699, 1473 

Kawaler, S. D. 2004, in Stellar Rotation, Proc. IAU Symp., 215, 

ed. A. Maeder k P. Eenes, 561 
Melis, C, Jura, M., Albert, L., Klein, B., k Zuckerman, B. 2010, 

ApJ, 722, 1078 
Rafikov, R. R. 2011a, ApJ, 732, L3 ( Rll) 
Rafikov, R. R. 2011b, arXiv:1102.4343 
Shakura, N. I. k Sunyaev, R. A. 1973, A&A, 24, 337 
Zuckerman, B. k Becklin, E. E. 1987, Nature, 330, 138 
Zuckerman, B., Koester, D., Melis, C, Hansen, B. M. S., k Jura, 

M. 2007, ApJ, 671, 872 
Walker, M. A. k Meszaros, P. 1989, ApJ, 346, 844 



APPENDIX 



FORMATION OF SHARP OUTER EDGE. 



Equation ([20]) is similar to Hopf (or simple wave) equation dy/dT + ydy/dx — 0, which is known to result in 
appearance of multivalued solutions usually interpreted as evidence for the shock formation. We may expect something 
similar in our case as well because the nonlinearity is clearly present in equation (|20[) . The equation for characteristic, 
crossing the point (x±, y = yi), is 



xy + ln(l — xy) = x\y + ln(l — x\y) + y 2 T. 
Since the left hand side of (|A1[) is negative for all xy, characteristic does not exist at times 

x\y + ln(l — X\y) o , , , it — 1 



T>-- 



(Al) 



(A2) 



y ~ (1 - e-") 2 

From that we can derive an upper limit on the time at which formation of a sharp edge occurs: 

T br eak < min[x 2 cp(r\\ (x))] Xg/2 ~ T thln < T thlck . (A3) 

Thus, we expect a discontinuity in surface density to form during disk evolution and our numerical results in ij44.ll 
support this expectation. 



